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Abstract 

A solid object can be coated by a nonwetting liquid since a receding contact line cannot exceed a 
critical speed. We theoretically investigate this forced wetting transition for axisymmetric menisci 
on fibers of varying radii. First, we use a matched asymptotic expansion and derive the maximum 
speed of dewetting. For all radii we find the maximum speed occurs at vanishing apparent contact 
angle. To further investigate the transition we numerically determine the bifurcation diagram for 
steady menisci. It is found that the meniscus profiles on thick fibers are smooth, even when there is 
a film deposited between the bath and the contact line, while profiles on thin fibers exhibit strong 
oscillations. We discuss how this could lead to different experimental scenarios of film deposition. 
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I. INTRODUCTION 



A convenient way to deposit a thin liquid layer on a surface is by withdrawing a solid 
from a liquid reservoir. The film is dragged along with the solid due to the viscous friction 
of the liquid. This principle is known as dip-coating and is a commonly used technique in 
industrial contexts l|, |2| . Once deposited on the surface, the film often has a thickness as 
predicted by Landau, Levich jsl and Derjaguin j4|, scaling with speed U of withdrawal as 
h oc f/^/^. Recently, however, a different class of solutions were identified, which are much 
thicker and scale as h cc U^^"^ jsl- These thick films were indeed realized experimentally in 
the case where the solid was partially wetting. 

The conditions of partial wetting introduces another interesting feature, namely that 
the film entrainment only appears above a critical velocity of withdrawal |6l-ll0| . Below 
this speed the contact line finds at a steady position, indicated as the meniscus rise A 
(Fig. H]). Due to viscous drag between the liquid and the solid, the dynamical position of 
A is higher than at equilibrium. This means that the apparent contact angle 9ap of the 
dynamical meniscus is smaller than the equilibrium angle Of.- The simplest interpretation of 
the transition to film deposition is that the apparent contact angle 9ap — )■ at some critical 
plate velocty. This idea was already postulated by Derjaguin and Levi fl. although the 
energy argument given by de Gennes |q] suggested a nonzero Oap at the transition. The 
hypothesis of ^ap=0, however, was given a rigorous mathematical basis (for a flat solid) by 
asymptotic expansions of the lubrication equations 3]. Actually, it was shown by 10 1 
that de Gennes energy argument can be extended to incorporate interface curvature: this 
exactly gives the lubrication equation, meaning that also the energy argument leads to a 
zero Oap at the transition. This theory gives a simple prediction for the maximum rise, based 
on the static meniscus solution with vanishing contact angle - for a fiber of radius tq this 



simply becomes 
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15| 



ro(ln^-c) forro<4 
a/24 for ro > 

Here ^c = {l/pgY^"^ is the capillary length based on surface tension 7, density p, gravity g 
and c is Euler's constant (0.57721). At intermediate radii tq ~ 4, the maximum rise can be 
determined numerically. 

2 




Uo 



To J 




X 


A 






FIG. 1: Schematic representation of the dip-coating setup: A fiber or cylinder of radius ro is 
withdrawn with speed Uq from a bath of viscous hquid. The axisymmetric meniscus profile is 
characterized by h{x), while A is denotes the maximum rise above the reservoir. 

Experimentally, the description of the forced wetting transition has remained ambiguous. 
The condition of a vanishing apparent contact angle was convincingly shown by Sedev & 
Petrov 16|. When withdrawing fibers or thin cylinders (ro/i'c ~ 0.06 — 1), they found a 
maximum rise of the meniscus consistent with ([T]) . Using cylinders of larger radii (ro /^c ~ 
10), Maleki et al. 17| found zero or nonzero 9ap at the transition, depending on the way 
9ap was determined. When using the criterion based on the meniscus height, the transition 
was found slightly before reaching Amax- Yet another set of experiments using a fiat plate 
(^n/^r, = oo) displayed a transition to film deposition clearly before reaching the maximum 



rise 
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19| . Still, during the unsteady entrainment phase the maximum recorded speed was 
reached exactly at v^^c- Note that in these experiments, the deposited liquid was not simply 
the Landau-Levich-Derjaguin film, but gave rise to thick films and even shock solutions. It 
was argued that the presence of these dynamical solutions are related to the pre-critical 



onset of entrainment 20|], but an explanation is still lacking. 

An additional complexity is that the contact line can spontaneously develop sharp corner 
structures, or even zig- zags. This has been observed in dip-coating 2lJ, splashing 22| . 



immersion lithography 23|, \2S\ and for drops sliding down an inclined plane 25, 

.Q- The 

conical structure of the interface near the contact line renders the problem truly three- 



dimensional, which affects the balance of the capillary forces 



27j . For sliding drops, it has 



been observed experimentally and described by a 3D lubrication model, that this change 
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in geometr y in deed leads to a nonzero apparent contact angle at the transition to liquid 



deposition 
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291]. This raises the question of how the geometry of the flow can influence 



the critical speed of wetting [30|. 

In this paper we theoretically study the withdrawal of fibers of arbitrary radii. By varying 



the ratio ro/^c, we continuously cover the range from thin fibers to the fiat p 
extend the asymptotic analysis that was previously done for the fiat plate 



ate. First, we 
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13| to the 



limit o: 
line 31 



thin fibers (Sec. |TT])- To resolve the singularity of viscous stress near the contact 



Si, we introduce^a slip length XUM- Other types of .mcroscopic regularization 
will give similar results [10']. Typical values for the slip and capillary lengths are A ~ 10 m 



and ic ~ 



10 ^m respectively. We can thus exploit the hierarchy of length scales 



A < ro < ic, (2) 

and perform a matched asymptotic expansion. The control parameter is the capillary number 
Ca = UoTj/j, which is the speed of withdrawal scaled by viscosity t] and surface tension 
7. The analysis yields the critical capillary number, which depends on the value of tq, 
and confirms that the maximum speed coincides with 6ap = 0, for all fiber radii tq. In 
this sense, the change in geometry does not qualitatively change the nature of the critical 
point. However, striking differences do show up when computing numerically the complete 
bifurcation diagrams for all steady solutions (Sec. IIII|) . These diagrams include solution 
branches above Amax that are unstable, but which have been observed as transients during 
film deposition for the plate case 191]. We find that for small fiber radii much below ic, the 
steady solutions no longer smoothly join the film solutions that mediate the deposition. In 
the Discussion section we speculate that this is why, experimentally, it is easier to approach 
the critical point for thin fibers (Sec. II Vp . 



II. ASYMPTOTIC ANALYSIS 



We compute the shape of an axisymmetric meniscus on a fiber of radius tq using the 
method of matched asymptotic expansions. The interface is characterized by h{x), as 
sketched in Fig. [H The matching procedure is outlined schematically in Fig. [2l At small 
scales, the dominant balance is between viscosity rj and surface tension 7, and is character- 
ized by the capillary number Ca. Viscous effects can be neglected on large scales, for which 
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FIG. 2: Schematic diagram showing the different asymptotic regions for the case of a thin fiber. The 
inner region originates from a balance between viscosity and surface tension. It has a microscopic 
contact angle 9(.. The outer region is a static meniscus joining a fiber with an apparent contact 
angle Oap- When the fiber radius tq ^ ^c, the outer profile is further separated into two regions [l^. 

the interface profile is that of a static meniscus. The problem is closed by matching the 
inner and outer solutions. The analysis provides the meniscus rise A as a function of Ca as 
well as the critical speed, both of which can be observed experimentally. We consider both 
large fiber radii (ro ^ d-c) and small fiber radii {tq <^ ic)- In all cases we take ro and ic to 
be macroscopic and much greater than the microscopic cutoff. Throughout the analysis, we 
scale all lengths by the capillary length, i.e. £c = 1- 



A. Inner solution: lubrication approximation 

To distinguish the solution h in the inner region and the outer region, we denote hin{x) 
as the solution in inner region and hout{x) as the solution in outer region. The characteristic 
length scale for the inner solution comes from the cutoff of the viscous singularity, which 
here we take the slip length A. Since typical interface curvatures turn out ~ Ca}^'^/X, as 
can be observed from the rescalings below, we can neglect the curvature contribution due to 
axisymmetry, which is of order I/tq. Hence, for the inner solution we can follow the analysis 



by Eggers 12|, |l3|, which was originally derived for the flat plate. For completeness, we 
briefly summarize the analysis and the central results. 

By restricting the analysis to small contact angles, h^^{0) = 9e I, one can determine 



h{x) from the lubrication approximation 35] : 
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Since the slip length A is the only length scale, we rescale the solutions according to 



h.4^)-3XH{0), (4) 

Hence @ reduces to 

where we introduced a reduced capillary number 6 = 3Ca/6'g. The boundary conditions are 

Hi^ = 0) = 0, (6) 

H'i^ = 0) = 1 (7) 

and the asymptotic behavior that has to be matched to the outer solution. Away from the 
contact line, where H ^ 1, further reduces to 

y'" = (8) 

y 

where we have put = ^^^^viO- This equation has an exact solution, whose properties 



have been summarized in 



361]. In parametric form, a solution with y{0) = reads 



l3{aAi{s)+l3Bi(s)) 
1 



s e [si,oo[, (9) 



y ^ (aAi(s)+/3Bi(s))2 

where Ai and Bi are Airy functions \sT\. The limit ^ — )■ corresponds to s — )■ oo, the 
opposite limit ^ — t- oo to s — )■ Si, where Si is a root of the denominator of ([9]): 

aAi{si) + PBi{si) = 0. (10) 



Since the solution extends to s = oo, Si has to be the largest root of ffTOj) . 

The solution ?/(^) is thus characterized by three parameters a, (3 and Si. Note that these 
are related according to ffTOj) . so that only two parameters are independent. The constant /3 
can be determined by matching (Q, which is valid only for ,^ > 1, to a solution of ([5]), which 



includes the effect of tfie cutoff and is tlius valid down to the position ^ = of the contact 
line [13|. It was found that 



13^ = 7rexp(-l/(35))/22/3 + 0(5), (11) 

which eliminates one of the two free parameters. The remaining parameter will be eliminated 
below by matching the large scale asymptotics of y{^) the outer solution of the problem. 
For that, we only need the asymptotic behavior of y{E,) for large ^, which reads: 

y{i) = \^ye + hyi+0{l), (12) 

where 

B. Outer solution: static meniscus 

At the scale of outer solution one can neglect viscous effects, and the profile is governed 
by surface tension and gravity. Thus equating the hydrostatic pressure and the capillary 
pressure gives 

K = A - X, (14) 

where k is the curvature of the interface. Remind that we expressed all lengths in the 
capillary length = 1- The curvature can be expressed from the geometric relation 

h" 1 

(l + Ct)^/^ (ro + /^out)(l + 0^/'' ^ ' 

The corresponding outer solution hout{x) is that of a meniscus of a liquid reservoir joining 

the fiber surface. The contact angle of the meniscus at the surface is denoted as the apparent 

contact angle, 6ap, since it refers to the apparent angle on the scale of the outer solution. 

The boundary conditions therefore are: 





= 0) 


= 0, 


^out (-^ 


= 0) 


a 

"api 


^out(-^ " 


= A) 


= OO. 



(16) 
(17) 
(18) 



For the present analysis we require only the asymptotic behavior near the contact line, which 
is obtained by a Taylor expansion, 

Kut{x) =eapX+^KapX^ + 0(^X^) . (19) 

Note that we consider small 6ap, since the inner solution is obtained in the lubrication limit. 

In general, the governing equation (1141) cannot be solved analytically. In the following we 
will consider two extreme cases for which analytical solution can be obtained, namely the 
larger fiber radius case (ro ^ 1) and the small fiber radius case (ro <^ 1). 



1. Large fiber radius: ro ^ 1 

In the case where the fiber radius is much larger than the capillary length, the second 
term on the right hand side of f[T^ due to the curvature of the fiber can be neglected. Then 
f lT^ can be written as 

- A - X. (20) 



(1 + hi,r/^ 

Integrating fl20|) once with respect to a;, we obtain 

where the boundary condition /i^^^ — t- oo at the position of the reservoir (x = A) is used. 
Evaluating ( l2Tl) at the contact line position [x = 0) and using the geometrical connection 



sin 9 = /iout/ V 1 + ^out) we end up with 



and f lT4|) immediately gives 



A = ^2(1- sin M (22) 
~ V2{1 - eap/2) (23) 

Kap-V2{l-eap/2). (24) 



2. Small fiber radius: ro <C 1 



For thin fibers it has been shown that the outer region can be further divided into two 
subregions 15|, as has been sketched in Fig. |2j In the region far away from the fiber [h^ 1), 
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the term due to the curvature of the fiber can be neglected. On the other hand, gravity can 
be neglected in the region close to the fiber {h '^1), and the meniscus is determined by the 
balance between the two curvature terms in fllSp . The profile near the fiber is a classical 
zero curvature interface that can be expressed as 



cosh 



) + sin Oap sinh f 

Jo cos 9ap J \rocos9, 



- 1 



ap f 



(25) 



In the following paragraphs we will match this small-scale part of the outer solution to the 
viscous inner solution. We therefore make a Taylor expansion, for small values of Qap^ 



2ro 



(26) 



To express our results in terms of the meniscus rise A, we quote the result obtained by 



James 



15| in which the two subregions of the outer meniscus were matched: 



A = ro 



In 



ro(l + sinOap) 



(27) 



where c is Euler's constant (0.57721...). 



C. Matching 



We are now in a position to perform the matching between inner and outer solutions. 
First, we write the inner solution in terms of the original variables. 



+ byO^X + O (1) 



(28) 



Once more, we separately discuss the limits of large and small fiber radii. 



1. Large fiber radius: tq ^» 1 

Comparing the inner solution ( 128|) to the outer solution (11911241) one finds the matching 
conditions 
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FIG. 3: (color online) Apparent contact angle 9ap versus Ca (A = 10"*^, = 0.1 radian) for large 
radius (rg = 1000) and small radius (rg = 0.01). Curves: result from asymptotic matching, solid 
curve: rg = 1000, dashed curve: vq = 0.01; Symbols: numerical result, squares: tq = 1000, circles: 
ro = 0.01. 




Oap = S'%e,, (29) 
= V26'/'^. (30) 

Adding these two conditions leads to an equation for Si clS Sb function of 6: 

2/9, 2^/W(.i) _ 2V6exp[-l/(3^)] 

Once Si is known, one can compute the apparent contact angle 

A typical result for 9ap as a function of Ca is shown in Fig. [3] (solid curve: ro = 1000). At 
vanishing speed, one recovers the equilibrium contact angle {6, = 0.1 in this example). The 
apparent contact angle decreases for increasing speed, and tends to 9ap = at a critical value 
Cac. The prediction from the matching compares very well to direct numerical solution of 
the problem, which will be discussed in the following section (solid squares). Of course, it 
is also possible to determine the critical speed directly from ( 13T1) . as shown in [3]. The 
critical value 6c is obtained when the Airy function takes its global maximum, Ai'{si) = 0, 
corresponding to s^ax = —1-088 ■ • ■. This gives a critical speed 
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^H2^/^3vr(i(w))2AjJ • ^^^^ 
Note that b = SCa/O^ and Ai{s max) =0.53566.... Physically, this corresponds to a vanishing 
apparent contact angle, as can be seen from (l32l) since Ai'{si)=0. Indeed, this confirms the 



conjecture by Derjaguin and Levi ll| that the maximum speed is attained when 6ap = 0. 



2. Small fiber radius: ro <C 1 

We now perform a similar analysis for thin fibers by using the outer solution (^E^, which 
was never worked out previously. Comparing this to the inner solution f l28p one finds the 
matching conditions 



Gap = S^^^byOe, 



(34) 



1 

To 



3A ■ 



(35) 



The parameter si can be solved as function of 5 from fl5Sl) . More explicitly, we can write 
dSS} as 



22/3roexp[-l/(35)] 



(36) 



51/3 37rAi2(si)A/^e ■ 
The apparent contact angle follows from flMl) . Since this condition is the same for both 
small fiber radius and large fiber radius, the explicit form of 9ap is also given by eqn. (132|) . 

Once again, solutions of the matching conditions cease to exist at a critical speed, which 
occurs when the Airy function takes it global maximum, Ai'{si) = 0. In perfect analogy to 
the flat plate case, this corresponds to 9ap = 0. The critical speed is given by 



^''[2y^37r{M{sL.)y\)\ ■ ^^^^ 
This result has the same structure as (1331) . valid for tq ^ 1. Apart from numerical coef- 
ficients, the main difference is that the fiber radius tq appears inside the logarithm as the 
relevant outer length scale; for the flat plate the outer scale is the capillary length. 

This result is further illustrated in Fig. |3] showing 9ap for a radius tq = 10"^ (dashed 
curve). The curve is similar to that obtained for a plate of infinite radius, with a vanishing 
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contact angle at the critical point. Note that this critical speed depends weakly (logarith- 
mically) on the fiber radius, in agreement with prediction f l37|) . In addition, there is a also 
a logarithmic dependence of 5c on the equilibrium contact angle 9^. Let us emphasize that 
the validity of the asymptotic analysis requires \/9(. to be small. This means that, strictly 
speaking, we cannot deal with extremely small values of Be- 



lt is instructive to compare our results with Voinov's formula {g]. The prediction by 
Voinov for 5c has the same structure as ours, but the factor inside the logarithm is not 
precisely specified (a ratio between the macroscopic length scale and the microscopic scale). 
In fact the factor refiects the dependence on the specific geometry of the problem, that in our 
approach is determined by the matching of the inner region and the outer region. Naturally, 
the inner scale turns out to be the slip length, while the outer scale is the fiber radius or 
the capillary length. However, Voinov's formula misses details like the factors 6'e, 5c inside 
the logarithm. Also the resulting 9ap vs Ca is a bit different from Voinov's formula, as was 



previously discussed in detail by Eggers |l3 |. 



III. NUMERICAL SOLUTION 

We now perform a numerical analysis of the fiber withdrawal problem. This will confirm 
the validity of the asymptotics and extend the results to ro ~ 1. However, the main added 
value is that the numerical solution can determine the complete bifurcation diagrams of 
dewetting for arbitrary tq. These contain steady state solutions above A^ax that serve as 



transients towards film deposition 
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and thus provide crucial additional information. 
Below we first develop a lubrication model that accounts for the axisymmetric nature of the 
fiow. This quantitative correction with respect to the fiat plate will turn out important for 



the bifurcation diagram. We then summarize the numerical results. 
A. Lubrication approximation on a fiber 

To formulate a hydrodynamic model for the axisymmetric meniscus on a fiber, we consider 
Stokes equations 

- Vp + ■qV'^U - V<I> = 0, (38) 
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V ■ f7 = 0, (39) 

where p is the pressure field in the hquid, rj is the viscosity of the hquid, U is the velocity field 
in the frame comoving with the fiber, and $ is the gravitational potential per unit volume 
in the liquid. Since the meniscus is axisymmetric, the velocity in azimuthal direction is zero. 
We consider small contact angle, 6'e ^ 1, thus the fiow is mainly in the vertical x direction, 
namely, the radial component of velocity is much smaller than the vertical component (i.e. 
\Ur\ <^ \Ux\)- The fiow is solved with a no-stress condition at the interface, is located at 
r = tq + h, and reads (in the frame of the fiber) 



H^ + ^j -,^Uo.. = 0. (40) 

\ / r=ro+ri 

At the fiber surface, r = ro, we apply a Navier slip boundary condition 



TT I - (A1\ 
^x\r=ro Q \r=ro- \ / 



dr 

The axial (vertical) component of the velocity field then becomes 



X 



2 Vro 



1 9(p + $) 
21] dx 



\\2h + 



ro 



(42) 



For thin films h/vQ ^ 1 this reduces to the usual parabolic Poiseuille profile, but quantitative 
corrections appear when h/vQ ~ 1. 

The lubrication equation is obtained by imposing a zero fiux condition in the frame of 
the reservoir 



/ (f/x + Uo)rdr = 0. 

Jro 



(43) 



With this, ( l42l) can then be simplified as 



d{p + $) 
dx 

where we introduced d=h/ro and 



SriUofid) 



h[h + 3X{l + d/2)f{d)]' 



(44) 
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^ 8^' (2 + d) 

^ 3[4(l + rf)4 1n(l + d) -rf(2 + d)(2 + 6d + 3d2)]- ^ 

This function is a correction factor with respect to the flat plate {d = 0), and has the 

property /(O) = 1. Finally, we replace the pressure by the Young-Laplace equation, 

p-Po = -in, (46) 

where k is the curvature of the interface given by (fT5l) . This gives the lubrication equation 
on a fiber 



^ _ 3Ca/(rf) 

/i[/i + 3A(l + c//2)/((i)] ■ ^ ^ 

Note that once again all lengths are scaled by the capillary length. For d = h/ro <^ 1, we 
recover the usual lubrication equation since /(O) = 1. 



B. Results 

1. Critical speed 

The above lubrication equation (H7!) is solved numerically with boundary conditions 



h{0) = 0, (48) 
h'iO) = e,, (49) 

imposed at the contact line and 



h'{A) = oo, (50) 
k{A) = 0, (51) 

at the reservoir. We varied ro and A and determined the meniscus as a function of Ca. 

Figure m shows the meniscus rise A as a function of Ca on a fiber of radius ro = 10~^. 
Different symbols correspond to different values of the slip length. In all cases we find a 
critical Cac above which solutions cease to exist. This indeed occurs close to Amax corre- 
sponding to a vanishing 9ap, which is indicated by the horizontal dotted line. The curves 
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FIG. 4: (color online) Height of meniscus A versus speed for different slip lengths (ro = 10~^, 
6e=0.1 radian). Curves are results from asymptotic matching, solid curve (the most left curve): 
A = 10~*, dashed curve: A = 10^^, dotted curve (the most right curve): A = 10^^. Symbols are the 
corresponding numerical results. As the slip length gets smaller, the agreement between numerics 
and asymptotic matching becomes better. The horizontal dotted line indicates the maximum height 
of meniscus A = 0.05414 calculated by ([27]) with 6ap = 0. 

are the predictions from the matched asymptotics, showing a good agreement with the nu- 
merical solutions. In particular, one observes convergence as the slip length is reduced from 
A = 10^®, 10^^ to 10^*^. This is because the separation of the two length scales A and tq 
is enhanced, which improves the validity of the matching asymptotic expansion. The same 
results were previously reported in Fig. [3], expressed in terms of 9ap rather than A. 

It is interesting to show how the critical speed Ca.c depends on the fiber radius tq. The 
numerical results are plotted as squares in Fig. O In agreement with the asymptotic analysis 
one observes two regimes. At small radii, ro ^ 1, the critical speed depends logarithmically 
on the radius. The solid red line is the asymptotic result fl37j) . For large radii the speed 
approaches the value of the fiat plate (133|) . indicated as dashed black hne. Indeed, the 
cross-over occurs for fibers with a radius that is comparable to the capillary length tq ~ 1. 

2. Meniscus rise: bifurcation diagram 

The results shown in Fig. H] represent only the lowest branch of solutions of a more com- 
plete bifurcation diagram. Indeed, one can identify solutions with A extending to arbitrary 
height above the meniscus, which are all characterized by Ca < Cac. These are summarized 
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FIG. 5: (color online) Critical capillary number Cac versus fiber radius ro (^e=0.1 radian, A=10^^). 
Squares: numerical results; Curve: Result from asymptotic matching for small fiber radius [Eq. 
()37p ]: Dotted line: Result from asymptotic matching for large fiber radius [Eq. (|33p ]. 



in Figs, in^ce for different fiber radii (all curves correspond to 6'e=0.05 and A=10~^). For 
ro = 1000 (Fig. E^l), we see after reaching Cac, the curve turns back to Ca < Cac but 
with solutions of increasing A. We refer to these solutions as the second branch, which is 
known to be unstable 20|. Further upwards we observe a series of bifurcations to higher 



branches, oscillating around a characteristic value Ca*. Typical meniscus profiles are shown 
in Fig. [6)d - in order to compare the profiles we have shifted the positions of the contact 
line such that the baths collapse. Following the bifurcation diagram, the profiles evolve to a 
film solution for which the contact line has moved to arbitrary height above the meniscus, 
i.e. A — )■ oo. This film solution (shown in Fig. Od) is not the Landau-Levich-Derjaguin film, 
but corresponds to the new class of "thick film" solutions identified in 

We then decrease the fiber radius to tq = 0.1, as shown in Figs. [Htd. We find that 
Cac increases almost by a factor 2 with respect to the large radius. By contrast, Ca* 
corresponding to the thick film increases only by a small amount. As a result the values 
of Cac and Ca* have become more separated. Also, the corresponding meniscus profiles 
display more structure. The thick film exhibits much stronger oscillations before joining the 
reservoir. These trends becomes more dramatic up further decreasing the radius ro = 0.001 
(Figs. [6^d). The difference between Cac and Ca* is very pronounced, and A changes much 
more dramatically for the 3rd branch solutions. In this sense, the bifurcation diagram has 
a very different structure from those of large fiber radius. Interestingly, there still exists a 
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FIG. 6: (color online) Bifurcation diagrams of steady solutions. Panels (a), (c), (e): bifurcation 
diagrams for ro=1000, 0.1 and 0.001 respectively. All curves correspond to 9(, = 0.05 and A=10^^. 
Panels (b), (d), (f): interface profiles for ro=1000, 0.1 and 0.001 respectively. We report profiles 
from the 2nd branch and 3rd branch, corresponding to solutions indicated in the bifurcation di- 
agram by large circles. The "thick film" solutions correspond to profiles "without contact line, or 
A —7- oo, for "which we define Ca = Ca*. The thickness of the film for tq = 0.001 is shown by the 
horizontal line just above x axis in (f). 

thick film solution matching to the bath, but the profile displays many oscillations (Fig. Ef). 
These oscillations decay only very slowly "when moving further away from the bath - the 
asymptotic thickness of the film is indicated by the horizontal line just above x axis. While 
for the flat plate the thick film solutions have been observed experimentally Isl, we expect 
the oscillatory solution obtained for small radii to be unstable and of no physical relevance. 

For completeness, we report the values of Cac and Ca* for different radii in a separate 
graph (Fig. [7j). Note that the theoretical curve for Cac deviates from the numerical results 
as early as vq ^ 10^^. The reason is that here we use a realistic value for the slip length 
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^0 

FIG. 7: Cac and Ca* as function of rg for the same parameters used in Fig. [6l Triangles: numerical 
results for Cac; Squares: numerical results for Ca*; Curve: Result for Cac from asymptotic matching 
for small fiber radius [Eq. (|37p ]: Dotted line: Result for Cac from asymptotic matching for large 
fiber radius [Eq. ([33]) ]. 

A = 10~^ (corresponding to ~ lOnm), instead of A = 10~^ used in Fig. [51 Clearly, the scale 
separation required for the asymptotic analysis starts to break down when the ratio A/ro is 
no longer very small. 



IV. DISCUSSION 



We investigated the steady-state profiles of axisymmetric menisci on a fiber that is with- 
drawn from a viscous liquid. The main motivation for this work was the mixed experimen- 
tal observations on the transition to film deposition obtained for fibers, large cylinders and 
plates. Sedev & Petrov [16^ found that the maximum steady profile has a meniscus rise iden- 
tical to a perfectly wetting liquid at equilibrium, suggesting a vanishing apparent contact 



ang. 



e O np- Other experiments found that steady-state solutions disappeared at a nonzero 
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n 



181], although the critical point could be accessed during transients [l9|]. Our present 
calculations show that steady solutions always cease to exist at Onp = 0, independent of 



the fiber radius. In addition, stability arguments put forward in 



13 



l38| suggest that all 



solutions of the lowest branch are perfectly stable up to the maximum speed, consistent with 



a saddle-node bifurcation 



20l ]. In that sense, our results do not provide an explanation why 



experimentally it is practically impossible to achieve steady menisci closer to 9ap = 0. The 
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main effect that was not talcen into account in our calculations is contact angle hysteresis 
due to heterogeneity of the substrate [sol, which was previously suggested to affects the 
details of the transition 40|. It has remained a challenge, however, to incorporate this into 
a full hydrodynamic description of moving contact lines. Essentially one has to modify the 
boundary condition by imposing a time-dependent microscopic contact angle at the moving 
contact line. 

The bifurcation diagrams calculated in the second part of the paper, however, do provide 



a new experimental perspective on the dynamics of film deposition. As shown in [19|], such 
bifurcation diagrams may be probed experimentally as transient states during entrainment. 
Namely, for Ca > Cac the evolution of the meniscus exactly follows the bifurcation diagram 
when plotting A versus the relative contact line velocity with respect to the solid. For very 
large fiber radius the profiles with large capillary rise are smoothly connected to the bath 
by a film that only displays a small "dimple" close to the bath. These dimple solutions 
have indeed been observed experimentally when plates are withdrawn with speeds above 
the critical speed. By contrast, for small fiber radii these solutions exhibit very strong 
oscillations (Fig. [6]) and we expect these solutions to be very unstable. In that case another 
dynamical mode must appear in order to deposit a liquid film - for example, one could think 
of the classical dewetting rim at the contact line connected to a Landau-Levich film j4lj]. A 
further investigation of these transients above the critical speed, in particular for different 
radii, should give a more complete picture of the forced wetting transition. 
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